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Summary 

We consider an infinite 3-dimensional elastic continuum whose material points 
experience no displacements, only rotations. This framework is a special case 
of the Cosserat theory of elasticity. Rotations of material points are described 
mathematically by attaching to each geometric point an orthonormal basis which 
gives a field of orthonormal bases called the coframe. As the dynamical variables 
(unknowns) of our theory we choose the coframe and a density. We write down 
the general dynamic variational functional for our rotational theory of elasticity, 
assuming our material to be physically linear but the kinematic model geometrically 
nonlinear. Allowing geometric nonlinearity is natural when dealing with rotations 
because rotations in dimension 3 are inherently nonlinear (rotations about different 
axes do not commute) and because there is no reason to exclude from our study large 
rotations such as full turns. The main result of the paper is an explicit construction 
of a class of time-dependent solutions which we call plane wave solutions; these are 
travelling waves of rotations. The existence of such explicit closed form solutions is a 
nontrivial fact given that our system of Euler-Lagrange equations is highly nonlinear. 
In the last section we consider a special case of our rotational theory of elasticity which 
in the stationary setting (harmonic time dependence and arbitrary dependence on 
spatial coordinates) turns out to be equivalent to a pair of massless Dirac equations. 



1. Introduction 

We work in 3-dimensional Euclidean space and view it as an elastic continuum whose 
material points can experience no displacements, only rotations, with rotations of different 
material points being totally independent. Rotations of material points of the 3-dimensional 
elastic continuum are described mathematically by attaching to each geometric point an 
orthonormal basis. This gives a field of orthonormal bases called the coframe. 

The purpose of our paper is to develop a theory of elasticity on rotations, i.e. a theory of 
elasticity in which the coframe plays the role of the dynamical variable (unknown quantity) . 
Recall that in classical elasticity the vector field of displacements is the dynamical variable. 

Our motivation for studying such a seemingly exotic problem comes from three main 
sources. 

The first source is Cosserat elasticity. In 1909 the Cosserat brothers proposed a theory of 
elasticity (1) which generalised classical elasticity by giving each material point rotational 
degrees of freedom. Cosserat elasticity has since become an accepted part of solid mechanics, 
though for most real life materials effects resulting from rotations of material points are 
small compared to effects resulting from displacements. From a purely mathematical point 
of view classical elasticity and rotational elasticity are two limit cases of Cosserat elasticity. 
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One of these limit cases, classical elasticity, has been extensively studied so it seems natural 
to examine now the other limit case. 

The second source is teleparallelism (= absolute parallelism = fcrnparallclismus), a 
subject promoted by A. Einstein and E. Cartan (2, 3, 4) in the late 1920s. The idea of 
rotating material points lies at the heart of teleparallelism and can easily be traced back to 
Cosserat elasticity: when in 1922 Cartan started developing what eventually became modern 
differential geometry he acknowledged (5) that he drew inspiration from the 'beautiful' work 
of the Cosserat brothers. The relationship between Cosserat elasticity and teleparallelism 
is examined in detail in the review paper (6). 

The third source is the theory of liquid crystals and, in particular, the concept of an 
Ericksen fluid. According to (7), in a liquid crystal one can observe 'orientation waves 
which propagate, inducing little or no motion of the fluid' and Ericksen's mathematical 
model is the natural way of describing this phenomenon. The only difference between 
Ericksen's model and ours is that in Ericksen's model one attaches to each geometric point 
a single unit vector rather than an orthonormal basis, as we do. 

Our paper has the following structure. In Section 2 we define our dynamical variables 
(unknowns of our theory), in Section 3 we write down the kinetic energy and Section 4 
we write down the potential energy. The Lagrangian of rotational elasticity is written 
down in Section 5. In Section 6 we reformulate our model in the language of spinors 
and in Section 7 we discuss the corresponding Euler-Lagrange equation. In Section 8 we 
construct an explicit class of solutions which we call plane wave solutions; this construction 
is summarised in Theorem 8.3 which is the main result of our paper. Finally, in Section 9 
we compare our model with the massless Dirac equation. 

2. Setting the playing field 

We work in Euclidean space R 3 equipped with Cartesian coordinates x a , a = 1,2,3, and 
standard Euclidean metric. We denote time by x°. Partial differentiation in x° and x a , 
a = 1,2,3, is denoted by d and d a respectively. 

The coframe d is a triple of orthonormal covector fields & , j = 1, 2, 3, in R 3 . Each covector 
field ffl can be written more explicitly as "& 3 a where the tensor index a = 1, 2, 3 enumerates 
the components. The orthonormality condition for the coframe can be represented as a 
single tensor identity 

g = 8 jk W®ti k (2.1) 

where S is the Kronecker delta and g = g a p = 8 a p is the Euclidean metric. For the sake of 
clarity we repeat formula (2.1) giving tensor indices explicitly and performing summation 
over frame indices explicitly: 5 a p — i? 1 ^ 1 /? + "d 2 a ^ 2 p +$ 3 a $ 3 i3 where a and (3 run through 
the values 1,2,3. We view the identity (2.1) as a kinematic constraint: the covector fields 
are chosen so that they satisfy (2.1), which leaves us with three real degrees of freedom 
at every point of M 3 . 

We work only with coframes which have positive orientation, i.e. which satisfy the 
condition 

dett^a = +1 > 0. (2.2) 

If one views & a as a 3 x 3 real matrix-function, then conditions (2.1) and (2.2) mean 
that this matrix-function is special orthogonal. Thus, the coframe can be thought of as a 
field of special orthogonal matrices. 
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As dynamical variables in our model we choose the coframe ■& and a positive density p. 
Our coframe and density are functions of Cartesian coordinates x a , a — 1, 2, 3, as well as of 
time a; . At a physical level, making the density p a dynamical variable means that we view 
our continuum more like a fluid rather than a solid. In other words, we allow the material 
to redistribute itself so that it finds its equilibrium density distribution. Observe that the 
total number of real dynamical degrees of freedom contained in the coframe i? and positive 
density p is four, exactly as in a two-component complex-valued spinor field. 

Note that there is nothing wrong in taking a prescribed density (as opposed to a density 
which is a dynamical variable): the theory one gets is very similar to the one described in 
the current paper and most formulae carry through with minimal changes. 

Below is the list of the main assumptions on which our model will be based. 

Assumption 1: our model is geometrically nonlinear. This means that we do not linearise 
rotations and we do not linearise the density. In other words, we allow our material points 
to experience full turns and we allow our density to experience changes comparable to the 
density itself. 

Assumption 2: our material is physically linear. This means that our potential energy is 
chosen to be quadratic in torsion (the latter serves as the measure of rotational deformations, 
see subsection 4.1). Note that physical linearity does not contradict geometric nonlincarity: 
locally (in space and time) material points "do not know" that they may eventually 
experience full rotations and the density "does not know" that it may eventually experience 
a change comparable to its current value. 

Assumption 3: our material is homogeneous and isotropic. Homogeneity means that 
physical properties of the material are the same at all points of our continuum and isotropy 
means that there are no preferred directions. 

Assumption 4-' our model is invariant under rigid rotations of the coframe. By a rigid 
rotation of the coframe we understand the transformation 



where k is a constant special orthogonal matrix. The thinking here is that when we 
attach an orthonormal basis to each geometric point of our continuum there is no reason 
to associate one particular direction with t? 1 , another with i) 2 and a third with z? 3 . What 
matters is how these directions change when we move from one point to another, i.e. how 
orthonormal bases at different points differ relative to each other. A rigid rotation of the 
coframe means that we simultaneously rotate all our orthonormal bases by the same angle 
around the same axis. We view rigid rotations of the coframe as gauge transformations and 
assume that our model does not feel them. See also (8) for a detailed exposition of gauge 
theory for problems similar to the ones considered in our paper. 

3. Kinetic energy 

Kinetic energy is given by the formula 



where c kln is some positive constant and w is the (pseudo)vector of angular velocity 



W i — y 3 k & 



(2.3) 




(3.1) 



(3.2) 
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Here A is the exterior product and * is the Hodge star (A.l). 

In writing the formula for kinetic energy (3.1) we think of each material point as a uniform 
ball possessing a moment of inertia and without a preferred axis of rotation. 

We give for reference a more explicit version of the formula for angular velocity (3.2): 

co a = - J2\ «o*'i - $ j id W 3 . (3.3) 

4. Potential energy 

4.1 Measuring rotational deformations 

In order to write down the formula for the potential energy we need to measure deformations 
caused by rotations of the coframe. More specifically, we need to measure deformations 
caused by the fact that at different points the coframe is oriented differently Obvious 
candidates for a measure of deformations are the rank two tensors 

K*:=d&, j = 1, 2, 3, (4.1) 

or, in more explicit form, K^ a p :— d a ^ p- The problem is that taken separately the three 
rank two tensors (4.1) are not invariant under rigid rotations of the coframe (2.3). The 
natural way of forming a truly invariant object is to make one rank three tensor out of the 
three rank two tensors K ■? according to the formula 

K := SjkW <g> K k = SjkW ® dd k . (4.2) 

The rank three tensor K is invariant under rigid rotations of the coframe (2.3) and, 
moreover, the individual rank two tensors can be recovered from K as 1 g = ■&^ a K ai g 
so there is no loss of information. 

Let us examine the symmetries of the tensor K . Observe that formula (2.1) implies 

= d a9P7 = d a (5 jk ti j p'd k 1 ) = 5 jk {d a & l3 )$ k 1 + 5 3k W p (d a i3 k 7 ) = K iafj + K Pari 

which means that the rank three tensor K is antisymmetric in the first and third indices, 

K ia p = -Kp ai . (4.3) 

Now, let us introduce another rank three tensor 

T := SjkP ® dd k (4.4) 

where d stands for the exterior derivative. The tensor (4.4) is obviously antisymmetric in 
the second and third indices 

T a (j 1 = -T Q7/ 3 (4.5) 
and is expressed via our original deformation tensor (4.2) as 

T Q(37 = K at3l - K ai p. (4.6) 



Formulae (4.6) and (4.3) imply 
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where the last two identities were obtained from the first one by a cyclic relabelling of tensor 
indices. Adding up the first and second identities and subtracting the third one wc get 

K a fjy = -(T Q( 3 7 + T 7Ct( g — Tg 7Q ) = - (T a ^ + Tp aj + T ia p) (4.7) 

(here we also used (4.5)). Note that the argument carried out above is a rephrasing of 
the standard argument that for a metric compatible affine connection contortion can be 
expressed via torsion, see subsection 7.2.6 in (9). 

Formulae (4.6) and (4.7) show that the tensors K and T are expressed via each other 
so either of them can be used as a measure of rotational deformations. We choose to use 
the tensor T because it has a clear geometric meaning: it is the torsion of the teleparallcl 
connection generated by the coframe i9, see Appendix A of (10) for a concise exposition. 
An additional advantage of using the tensor T is that the definition (4.4) of this tensor 
does not require the use of covariant derivatives so it works when the metric g appearing in 
formula (2.1) is not assumed to be Euclidean. The latter was important for Einstein and 
Cartan who arrived at the mathematical model similar to the one described in in our paper 
coming from general relativity. Recall that in general relativity the metric plays the role 
of dynamical variable so for someone with a relativistic background assuming the metric to 
be Euclidean (i.e. space to be flat) is unnatural. 

Starting from Einstein's works (4) torsion is traditionally used as a measure of 
deformations when modelling elastic continua with rotations. We shall follow this tradition 
and construct our potential energy as a function(al) of T. However, before writing down the 
formula for potential energy we will simplify matters by using the fact that we are working 
in 3D (our previous arguments were dimension- independent). 

Applying the Hodge star (A.l) in the second and third indices we switch from the original 
torsion tensor T to the tensor 

T a p := ^T a i s e jS 0. (4.8) 

Of course, the tensor T can be recovered from T as 

T a /3j = T a espy. (4.9) 
* 

Formulae (4.8) and (4.9) show that the tensors T and T are expressed via each other so 
cither of them can be used as a measure of rotational deformations. We choose to use the 
tensor T because it has lower rank, two instead of three. 

Formulae (4.4) and (4.8) imply 



T = 8 jk d j ® *dd k = SjkP ® curltf' . 



(4.10) 
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We see that T is a rank two tensor without any symmetries and with arbitrary trace. This 
is the tensor we will be using as a measure of rotational deformations when writing down 

the formula for potential energy. The tensor T is sometimes called the dislocation density 
tensor (11) . 

We give for reference a more explicit version of formula (4.10): 

3 1021^3 -&> 1 d3 l d i 2 d j ld3&> 1 -'d j id 1 '& j 3 ■& j 1 d 1 ft> 2 -t^'iQ^'A 

T a0 = $ j 2d2$ j 3 - & j 2d 3 ti j 2 1^203^1 ~ ^2<9l^3 ^201^2 ~ 1^2 02^1 • (4.11) 

j=l ^302^3 - $ j 3d 3 $ j 2 ^383^1 - $ j 3dl$ j 3 & 2 - ^'s^'l/ 

4.2 Irreducible decomposition of rotational deformations 

Recall the logic of classical linear elasticity (12): after identifying the measure of 

deformation one decomposes it into irreducible pieces. We follow this logic by decomposing 

* 

the tensor T into irreducible pieces. The construction presented below is similar to (12), 
the only difference being that instead of a symmetric rank two tensor, strain, we deal with 

a rank two tensor, T, without any symmetries. 

Decomposing the rank two tensor T into irreducible pieces means the following. We fix 
a point in R 3 and at this point look at all real rank two tensors P. Such tensors can be 
viewed as elements of a real 9-dimensional vector space V equipped with inner product 

(P,Q)v :=P a pQ aP (4.12) 

and corresponding norm 

||P||v = V(P7PW = yJPa^. (4-13) 

Let us now examine what happens when we rotate our Cartesian coordinate system x a , 
i.e. when we perform a linear change of coordinates preserving the metric g a p and 
orientation. The components of our tensors P a p change in a particular way under rotations 
of the coordinate system, so we get an action of the group SO(3) on the vector space 
V. Looking for irreducible pieces of torsion means identifying subspaces of V which are 
invariant under the action of the group SO(3), i.e. which map into themselves, and which 
do not contain smaller nontrivial invariant subspaces. 

In our case the invariant subspaces are obvious. These are 

• the 1-dimensional subspace of real rank two tensors proportional to the metric, 

• the 3-dimensional subspace of real antisymmetric rank two tensors and 

• the 5-dimensional subspace of real symmetric trace-free rank two tensors. 

These three subspaces are clearly irreducible and mutually orthogonal in the inner 
product (4.12). 

Our rank two tensor T can now be written as a sum of three irreducible pieces 

* * * * 

rp Jiax _j_ rpVCC _j_ rptCYl ^ 14) 
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where 



T7 

-^g a p, (4.15) 

1 (4.16) 



* * * 

J a/3 — 1 <*P- 1 a [i~ 1 c[i ~ ^ ^~9a0- (4-17) 

We label the irreducible pieces (4.15), (4.16) and (4.17) by the adjectives axial, vector and 
tensor respectively, which is terminology traditional in alternative theories of gravity (6). 

4.3 Formula for potential energy 

Following the logic of classical linear elasticity (12) we now write down the explicit formula 
for potential energy: 



P(x°) = J (c ax ||T ax ||t> + c V0C ||T vcc ||t> + c tcn \\T tcn \\ 2 v )pdx 1 dx 2 dx 3 



(4.18) 



where c ax , c voc and c ten are some nonncgative constants (elastic moduli), not all zero, and 
| • ||y is the norm (4.13). Comparing our formula (4.18) with formula (4.3) from (12) wc 
see a difference between classical and rotational elasticity: classical elasticity involves two 
clastic moduli whereas rotational elasticity involves three. The extra elastic modulus c vcc 

is needed because the tensor T which we use as measure of rotational deformations is not 
necessarily symmetric. 

Formula (4.18) is the one traditionally used in teleparallelism. This formula already 
appears in the original papers of Einstein (4), though for some reason' 1 ' Einstein did not 

include the axial term c^HT^H 2 . Subsequent authors always used three terms, see, for 
example, formula (26) in (6). 

4.4 Simplifying the formula for potential energy 
Let us introduce the (pseudo)scalar 

f-=T a a (4.19) 

and the vector 

v a := T^e Ma . (4.20) 
Formulae (4.13)-(4.17), (4.19) and (4.20) imply 

llHlv = \f 2 ' ( 4 - 21 ) 
\\T™\\ 2 V = \\H\ (4.22) 



t The reason could be that Einstein was primarily interested in providing a geometric interpretation of 
electromagnetism and might have felt that the axial term would not contribute to the electromagnetic field. 
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\\T^r v = \\T\\ 2 v-lf-l\H 2 . (4.23) 
Substituting (4.11) into (4.19) and (4.20) we get more explicit formulae for / and v. 

3 

/ = ^(^i<9 2 ^3 - j ida# i 2 + ^2d 3 ^i - &2di& 3 + •0 i a&i# i 2 - ^s^'i), (4.24) 

Va = J2\ & j 3 d rf j 3 ~ Wadset - Wldi&2 + $ j id 2 # j i . (4.25) 
Substituting formulae (4.19)-(4.21) into formula (4.18) we get 

//ax _ ton vcc _ ten * * \ 

( ^—f + ~ 2"^— ll"l| 2 + c tcn T a/3 T a P J pdx x dx^dx z . (4.26) 

The advantage of writing potential energy in the form (4.26) is that the geometric quantities 

* 

/, w and T appearing in this formula have relatively compact explicit representations (4.24), 
(4.25) and (4.11). 

5. Lagrangian of rotational elasticity 

We combine our potential energy (4.26) and kinetic energy (3.1) in forming the action 
(variational function) of dynamic rotational elasticity 

S(0,p) = J (P(x°) - K(x°))dx° = J L(d,p)dx°dx 1 dx 2 dx 3 (5.1) 

where 

/ „ax _ ten vcc _ ten * * \ 

WP) = [ 3 f + 2 Ikll 2 + c tCD T a pT a P c kin \\oj\\ 2 j p (5.2) 

is the Lagrangian density. Recall that the geometric quantities f,v,T and to appearing in 
formula (5.2) are defined by formulae (4.24), (4.25), (4.11) and (3.3) respectively. 

Our construction of the action (5.1) out of potential and kinetic energies is Newtonian: 
compare with classical linear elasticity or even the harmonic oscillator in classical mechanics. 
An alternative approach is the relativistic one which boils down to rewriting the formula for 
potential energy in Lorentzian signature in dimension 1+3, with this "extended" potential 
energy becoming the action. The Newtonian and relativistic approaches are different which 
can be seen, for example, from the fact that the relativistic approach always imposes a 
unique velocity of wave propagation (speed of light) whereas with the Newtonian approach 
one expects to get at least two distinct wave velocities. 

Starting with Einstein, most authors working in the subject of teleparallelism adopt the 
relativistic approach. We shall, however, stick with the Newtonian approach (5.1). 
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6. Reformulating the problem in the language of spinors 

Our field equations (Euler-Lagrange equations) arc obtained by varying the action (5.1) 
with respect to the coframe ■& and density p. Varying with respect to the density p is easy: 

this gives the field equation ^-^ .f 2 + £— ^ \\v\\ 2 + c tcn T afj T af3 ~ c kin ||cj|| 2 = which 
is equivalent to L(i3, p) — 0. Varying with respect to the coframe d is more difficult because 
we have to maintain the kinematic constraint (2.1). 

This technical difficulty can be overcome by switching to a different dynamical variable. 
Namely, it is known (13) that in dimension 3 a coframe d and a (positive) density p arc 

equivalent to a 2-componcnt complex-valued spinor field £ = £ a = ^ 2 I modulo the sign 



e 

of £. The explicit formulae establishing this equivalence are 

p = eooa b £,\ (6.i) 

(tf 1 + id 2 ) a = p-VV o6a r<W d , (6.2) 

$\ = p-^tTaab?. (6-3) 

Here a are Pauli matrices and e is "metric spinor" (see (A.2)-(A.4)), the free tensor index 
a runs through the values 1, 2, 3, and the spinor summation indices run through the values 
1,2 or 1,2. The advantage of switching to a spinor field £ is that there are no kinematic 
constraints on its components, so the derivation of field equations becomes straightforward. 
We give for reference more explicit versions of formulae (6.1)-(6.3): 

P =ez i -ee. (6.4) 
fie? -(e? , 

(v 1 + ^v 2 ) a = (?e + eer 1 ^) 2 + *(a 2 1 , (6.5) 

V -zee 

( ~ee + 

$\ = (?e + eer 1 &e - 1 ■ (6-6) 




Let us rewrite the geometric quantities /, v, T and ui appearing in formula (5.2) in terms 
of the spinor field £. The spinor representation of angular velocity uj was derived in (13): 

. e« a abd a e - e<j a g b do^ 

w« = i To 77i ( 6J ) 



or, more explicitly, 



where the "c.c." stands for "complex conjugate term" . The spinor representation of the 
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tensor T is derived in Appendix B, see formula (B.l) or its more explicit version (B.7). 
Substituting (B.l) or (B.7) into (4.19) and (4.20) we arrive at spinor representations for the 
(pseudo)scalar / and vector v. 

f = - 2t ¥^ ' (6 ' 9) 

v a = -e^ (6.10) 

where d a := g al3 dp — d a , or, more explicitly, 



i 



i?d 3 e + i?d 3 e - i^o 1 ? + if d 1 ? | + c.c. (6.i2) 



Note that formula (6.9) is a rephrasing of formula (B.5) from (13). 

From now on we write our action (5.1) and Lagrangian density (5.2) as S(£) and 
rather than S($,p) and L($, p), thus indicating that we have switched to spinors. The 
explicit formula for L(£) is obtained by substituting formulae (6.1), (6.9), (6.10), (B.l) and 
(6.7) into (5.2). The nonvanishing spinor field £ is the new dynamical variable and it will 
be varied without any constraints. 



7. Euler Lagrange equation 

Let us perform a formal variation of our spinor field £ M> £ + <5£, where <5£ : R x R 3 — > C 2 
is an arbitrary (infinitely) smooth function with compact support. Then, after integration 
by parts, the variation of our action can be written as 

SS = J (FaS^ + F a 6£ a ) dx°dx 1 dx 2 dx 3 (7.1) 

where F is a dotted spinor field uniquely determined by the undotted spinor field £. The 
Euler-Lagrange for our unknown spinor field £ is, therefore, 

F = 0. (7.2) 

The map 

£ F (7.3) 

defines a nonlinear second order partial differential operator in the variables x° (time) and 
x a , a = 1,2,3 (Cartesian coordinates). 

We shall refrain from writing down the Euler-Lagrange equation (7.2) explicitly. The 
reason for this is that in the current paper we are interested in finding a particular class of 
solutions for which the procedure is much simpler. 
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Note that for the special case of a purely axial material, i.e. material with 

c voc = 0, c tcn = 0, (7.4) 

the Euler-Lagrange equation (7.2) was written down explicitly in (13). The calculations in 
(13) were carried out under the additional assumption 

kin _ ax 

" 3 

which can always be achieved by rescaling time x°. 

8. Plane wave solutions 

We seek solutions of the form 

£(x) = e" ipx C (8.1) 

where ( =/= is a constant (complex) spinor and p is a constant real covector. Here we 
use relativistic notation, incorporating time x° into our coordinates. This means that 
x = (x , x , x 2 , x 3 ) and p = (po,Pi,P2,Ps)', bold type indicates that we are working in 
(l+3)-dimensional spacetime. The number \pg\ is the wave frequency and the covector 
(pi,P2,P3) is the wave vector in original 3-dimensional Euclidean space. The 4-component 
covector p = (po,Pi,P2,P3) has the meaning of relativistic 4- momentum. 
Throughout this section as well as the next one we assume that 

Po + (8.2) 

which means that we are not interested in static (time-independent) solutions. The sign of 
Po can be arbitrary. 

Our Euler-Lagrange equation (7.2) is highly nonlinear so it is by no means obvious that 
one can seek solutions in the form of plane waves (8.1). Fortunately (and miraculously) this 
is the case. In order to sec this, we rewrite our Euler-Lagrange equation (7.2) in equivalent 
form 

Note that the sign in the exponent in (8.3) is opposite to that in (8.1). We have 

Lemma 8.1. If the spinor field £ is a plane wave (8.1) then the left-hand side of equation 
(8.3) is constant, i.e. it does not depend on x. 

Of course, Lemma 8.1 can be equivalently reformulated as follows: the nonlinear partial 
differential operator (7.3) maps a plane wave (8.1) into a plane wave with the same 
relativistic 4-momentum p. 

The proof of Lemma 8.1 is quite technical and is given in Appendix C. 

Lemma 8.1 justifies separation of variables, i.e. it reduces the study of the nonlinear partial 
differential equation (7.2) for the unknown spinor field £ to the study of the rational algebraic 
equation (8.3) for the unknown constant spinor £. We suspect that the underlying group- 
theoretic reason for our nonlinear partial differential equation (7.2) admitting separation 
of variables is the fact that our model is U(l)-invariant, i.e. it is invariant under the 
multiplication of the spinor field £ by a complex constant of modulus 1. Hence, it is feasible 



(7.5) 
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that one could prove Lemma 8.1, as well as Lemma 8.2 stated further down in this section, 
without performing the laborious calculations presented in Appendix C. 

We are now faced with the task of writing down the LHS of equation (8.3) explicitly and 
with minimal calculations. To this end we address a seemingly different issue: we examine 
what happens when we substitute our plane wave (8.1) into our Lagrangian density 
rather than the Euler-Lagrange equation (7.2). 

Substituting formula (8.1) into formulae (6.1), (6.9), (6.10), (B.l) and (6.7) we get 

P = ( a coab( b , 



/ = - 



T a (3 = 2 



= _2£^CV^o6PV 
K h O a abP0C h 



or, more explicitly, 
/ 4 



CY + C^C 2 , 



C 1 C 1 +C 2 C 2 



(M-c'e - c~V) + iP2(c\ 2 c~V) +M-c i c l + m) 



( 



* 

Ti 2 

T21 

^23 

* 

* 



C 1 C 1 +C 2 C 2 



P2 {c'e-c'ej+m (cv-cv 

\ipi (-c _i c 2 + cv) + P 2 (-c _i c 2 - cv 
( i P 2 (c _i c 2 - c 2 e) + P3 (-cY + c 2 c 2 ) ^ 

*Pi ( Jo • C~Y) 

pi f<Y-cY' 
^(cy + cY 

H 2 C 2 )+pi v 

p 2 fcY-cY 
/'.; (c~Y+cY 
(-c~Y+cY 



P3 



(-c~Y 



_ C 1 C 2 



cY) 



«P3 



-cY) 



(8.4) 
(8.5) 

(8.6) 

(8.7) 

(8.8) 

(8.9) 
(8.10) 

(8.11) 



(8.12) 
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Substituting formulae (8.4)-(8.8) or their more explicit versions (8.9)-(8.13) into formula 
(5.2) we arrive at a Lagrangian density L((;p) which does not depend on x. The self- 
contained formula for this Lagrangian density, written in terms of 4-momentum p and 
4-current 

ja:=C<wC 6 , (8-14) 

is 

L(C;p) = - (c voc +c tcn ) ||p|| 2 ||j|| 2 + - (V- V ec + \c'A (j. p f-±c ki »pl\\j\\ 2 (8.15) 
Jo Jo V 3 2 6 / Jo 

Here we write our 4-covectors as p = (po,p) and j = (jo, j), where p and j are 3-covectors. 

We view the 4-momcntum p as a parameter and the constant spinor £ ^ as the 
dynamical variable. The Lagrangian density L((; p) is a smooth function of Re£ and Im£, 
so varying £ we get 

5L = G a 5C + G a 5C (8.16) 

where G is a dotted constant spinor expressed via the partial derivatives of L(£; p) with 
respect to Re £ and Im£. It is natural to ask the question: what is the relationship between 
the spinor field F appearing in formula (7.1) and the constant spinor G appearing in formula 
(8.16)? The answer is given by 

Lemma 8.2. If the spinor field £ is a plane wave (8.1) then G = e ip x F. 
The proof of Lemma 8.2 is presented in Appendix C. 

Lemma 8.2 reduces the construction of plane wave solutions of rotational elasticity to 
finding the critical, with respect to £, points of the function L((;p). Varying (8.15), we 
arrive at the following equation for critical points: 

- (c vec + c tcn ) \\p\\ 2 j a a a a b C b - 4 ( cVCC + ctCn ) WpfWjfvoabt" 
Jo Jo 

+ 1 (|c- - \c- + \f) (J • P)P W - J (lc- - \c- + Jc*«) V • P fa 0ab C 

- ~c kin p 2 j a a a&b ( b + ^c kin po 2 \\j\\ 2 <Toab( b = 0. (8.17) 
Jo Jo 

Recall that the 4-current j = (jo,j) appearing in the above equation is defined in accordance 
with formula (8.14). 

It now remains to find the 4-momenta p and spinors C ^ which satisfy equation (8.17). 
We carry out the analysis of equation (8.17) assuming that 

P = Jo = CV + C 2 C 2 = 1. (848) 

Condition (8.18) is a normalisation of the density: general plane wave solutions differ from 
those satisfying condition (8.18) by a real scaling factor. Furthermore, we assume that 



(8.19) 
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Indeed, any spinor £ satisfying condition (8.18) can be transformed into the spinor (8.19) 
by the linear transformation 

C^UC (8.20) 

where U is a special (det U = 1) unitary matrix. The transformation (8.20) leads to a 
rotation of the spatial part of the 4-current (8.14), so plane wave solutions with general £ 
differ from those with ( of the form (8.19) by a rotation. 

Substituting (8.17), (8.19), (A.3), (A.4) and (8.13) into (8.16) we get 



4 (c voc + c ton ) b|| 2 (T) - 2 (c vec + c tol ) |bf (I 



0J > \0 

^ „vec i 1 „ten 1 „ / .P3 \ A ( ^ „ax 1 „vcc , 1 ten \ 2 f 1 



+ 8 o c ax - ^c vec + -c tcn p 3 , . - 4 -c ax - -c vcc + -c ten ^ . 
1 3 2 6 / y \Pi + IP2 J V 3 2 6 / V 



2 

kin / Po 



4C ( q ) = 0, 



or, equivalently, 



2 (c~ + c-) (11^+4 (§c- - ic- + Jc-) ( 2p3( ^ + . p2) )_ 4c - = 0. (8.21) 



Put 



4c ax + 2c tcn 



-,vec i «ten 



Wi: =V— U2:= V^r-- (8 - 22) 

Note that because we assumed our three elastic moduli to be nonnegative and not all zero, 
our vi and v 2 are nonnegative and not both zero. Using (8.22) we can now rewrite equation 
(8.21) in more compact form 

<^>«-« -(f)-* ^ 

The analysis of equation (8.23) is elementary and the outcome is summarised in the 
following theorem, which is the main result of our paper. 

Theorem 8.3. Plane wave solutions of rotational elasticity can, up to rescaling and rotation, 
be explicitly written down in the form (8.1), (8.19) with arbitrary nonzero po and p = 
(pi,P2,P3) determined as follows. 

• If vi > and w 2 > and v\ ^ u 2 then we have two possibilities: 

— p= ^0, 0, ±— ^ (type 1 wave), or 

— p = [ J^l cos m J^i s i n( p ; ) where ip G R is arbitrary (type 2 wave). 

• If vi > and v 2 > and «i = w 2 then p is an arbitrary 3- vector satisfying ||p|| = . 

• Ifu! > and w 2 = 0thcnp = (0,0, ±— J. 



II II 

• If v\ = and w 2 > then p = ( — — cos ip, — sin ip, I where ip 6 R is arbitrary. 

V 2 V2 
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Theorem 8.3 shows that rotational elasticity, like classical linear elasticity, produces two 
distinct types of plane wave solutions. We call these solutions type 1 and type 2 and they 
propagate with velocities Vi and v 2 respectively, with v\ and v 2 given by formulae (8.22). 

However, unlike with classical linear elasticity, in rotational elasticity the two wave 
velocities, v\ and v 2 , are not ordered, i.e. we do not know a priori which one, v\ or v 2 , 
is bigger. The reason the two wave velocities are not ordered is because rotational elasticity 
has three elastic moduli compared to the two elastic moduli of classical linear elasticity. The 
"extra" elastic modulus is c vcc , the one associated with the antisymmetric part of the rank 

two tensor T which we use as a measure of rotational deformations. If we set c voc = 0, we 

end up with the inequality v\ > \J\v 2 similar to the well known inequality from classical 

linear elasticity, see formula (22.5) in (12). 

9. The massless Dirac equation 

In this section we consider a purely axial material (7.4), (7.5). Substituting (7.4) and (7.5) 
into (8.22) we get v\ = 1 and v 2 — 0, so a purely axial material supports only type 1 waves. 
Throughout this section we also retain the assumption (8.2). 

Note that a purely axial material has the remarkable property that its potential energy 
is invariant under conformal rescalings of the spatial metric by an arbitrary positive scalar 
function. We do not elaborate on this issue in the current paper because we chose to work 
with a specific (standard Euclidean) metric. The appropriate arguments are presented in 
Section 2 of (13). 

Our aim is to compare our model with the linear partial differential equation (or, more 
precisely, system of two linear partial differential equations) 

i{±<y\ b d + <J a abd a )e = 0. (9.1) 

Here a are Pauli matrices (A. 3), (A. 4), the free spinor index a runs through the values 
1, 2, summation is carried out over the tensor index a = 1, 2, 3 as well as the spinor index 
b = 1,2, and £ is the unknown spinor field. We give for reference a more explicit version of 
equation (9.1): 

>{tS£ X -%){$)-"■ <-» 

Equation (9.1) is called the massless Dirac equation or the Weyl equation This equation 
is the accepted mathematical model for a massless neutrino field. The two choices of sign 
in (9.1) give two versions of the Weyl equation which differ by time reversal. Thus, we have 
a pair of Weyl equations. 

We want to compare plane wave solutions (see (8.1)) of our model with those of the Weyl 
equation. As both models are invariant under the rescaling of the spinor field by a positive 
real constant as well as the rotations of Euclidean 3-space, it is sufficient to compare plane 
wave solutions for ( of the form (8.19). Substituting (8.1) and (8.19) into (9.1) or its more 
explicit version (9.2) we get p = (0, 0, ±f>o) which is exactly what Theorem 8.3 gives us. 
Thus, we have established 

Theorem 9.1. In the case of a purely axial material a plane wave spinor field is a solution 
of rotational elasticity if and only if it is a solution of one of the two Weyl equations (9.1). 
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It turns out that, in fact, a much stronger result holds. Consider a spinor field of the 
form 

£(x°, x 1 , x 2 , x 3 ) = e- lpax ° ^{x 1 , x 2 , x 3 ). (9.3) 

We will call spinor fields of the form (9.3) stationary. In considering stationary spinor fields 
what we are doing is separating out time only as opposed to separating out all the variables. 
The following result generalises Theorem 9.1. 

Theorem 9.2. In the case of a purely axial material a nonvanishing stationary spinor field 
is a solution of rotational elasticity if and only if it is a solution of one of the two Weyl 
equations (9.1). 

Theorem 9.2 was proved in (13) and the proof is quite delicate. It involves an argument 
which reduces a nonlinear second order partial differential equation of a particular type 
to a pair of linear first order partial differential equations, which is, effectively, a form of 
integrability. An abstract self-contained version of this argument is given in Appendix B 
of (14). 
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APPENDIX A 

Notation 

Our notation follows (10, 13, 14, 15). The only difference with (10, 15) is that in the latter the 
spacetime metric has opposite signature. In (13, 14) the signature is the same as in the current 
paper, i.e. the 3-dimensional spatial metric has signature + + + . 

We use Greek letters for tensor (holonomic) indices and Latin letters for frame (anholonomic) 
indices. We identify differential forms with covariant antisymmetric tensors. 

We define the action of the Hodge star on a rank r antisymmetric tensor R as 

(*fl)« r+1 ...a 3 := (r!)- 1 fl ai - ap Ea 1 ...a 3 (A.l) 

where e is the totally antisymmetric quantity, £123 := +1. 

We use two-component complex-valued spinors (Weyl spinors) whose indices run through the 
values 1,2 or 1,2. Complex conjugation makes the undotted indices dotted and vice versa. 

We define the "metric spinor" 

ab dfr / 1 \ / A r»\ 

e ab = e dh = e =e = I 1 I (A.2) 



and choose Pauli matrices 

O~0a 



ab = = -0-°ab, (A.3) 

{0 A 1 /0 -A 2 A \ 3 , A A s 

Cld6 = I q I = O" ab, °~2ab = I ^ q I = CT ab, 0"3db = I Q _^ J = °~ ab ■ 

Here the first index enumerates rows and the second enumerates columns 

APPENDIX B 

Spinor representation of torsion 

We show in this appendix that the tensor T a p defined by formula (4.10) (it is the Hodge dual, 
in the last pair indices, of the torsion tensor) is expressed via the spinor field £ as 

Note that formula (B.l) is invariant under the rescaling of our spinor field by an arbitrary positive 
scalar function. 

Formula (B.l) is proved by direct substitution of formulae (6.2) and (6.3) into (4.10). In order 
to simplify calculations we observe that the expressions in the left- and right-hand sides of formula 
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(B.l) have an invariant nature, hence it is sufficient to prove formula (B.l) at a point at which the 

Then at this point we have 



spinor field takes the value 



where a = 1, 2, 3. Note that formulae (B.3) imply 



(B.2) 
(B.3) 



[curl^ 1 + M 2 )]p = 



-2d 2 e - daft 1 - if 1 ) 
29i£ 2 + d 3 (£ 1 -£ i ) 

M* 1 



(B.4) 



/ -d 3 (-ie+i?) 

[curl A, = ] (B.5) 

We now rewrite formula (4.10) in the form 

T = i^ 1 - M 2 ) <g) curl^ 1 + i?? 2 ) + ^(tf 1 + itf 2 ) <g) curl^ 1 - id 2 ) + tf 3 ® curli? 3 . (B.6) 



Substituting formulae (B.2), (B.4) and (B.5) into formula (B.6) we get 



V 



id 2 (e-e) 



which coincides with the RHS of formula (B.l). This completes the proof. 
We give for reference a more explicit version of formula (B.l): 



* 

Tl2 
Tl3 

* 
* 

T31 

* 



iL'hC - i^ 2 d 1 e 

i?d 2 e + i£}d 2 i 2 

ii ihe ■ <e<hi< 



+ c.c. 



(B.7) 
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APPENDIX C 

Separation of variables 

In this appendix we prove Lemmata 8.1 and 8.2. Note that it would suffice to prove Lemma 8.2 
only, because Lemma 8.1 follows from Lemma 8.2. However, we prove Lemma 8.1 first for the sake 
of clarity of exposition. 

Let us arrange the (pseudo)scalar /, the three components of the vector v a , the nine components 

* 

of the tensor T a p and the three components of the (pseudo)vector u) a into one 16-component 
"vector" Vj, J — 1, ... ,16. Then our Lagrangian density (5.2) can be written as 



L{i) = P Y J AjV] 



(C.l) 



where the Aj are some real constants. Put Wj := y/pVj. Then formula (C.l) takes the form 



(C.2) 



According to formulae (6.9), (6.10), (B.l), (6.7) and (6.1) the components of the "vector" W are 
expressed via the spinor field £ as 



Wj =i 



£, a B J gbda^ — £, b BJ ab d a £, a 

(£ £ <W) 1/a 



(C.3) 



where BJ a b are some constants and summation is carried out over the spinor indices a — 1,2, 
b = 1, 2, and over a = 0, 1, 2, 3. Here we use bold type to indicate relativistic notation, when time 
x° is viewed as one of the coordinates in (l+3)-dimensional spacetime. 

Note that for given J and a the 2x2 matrices BJ a b are Hermitian. This is because each of 
these matrices is a linear combination with real coefficients of the Pauli matrices a 13 ab , /3 = 1, 2, 3. 

Our action is S(£) = J L(£) where for the sake of brevity we dropped dx dx 1 dx dx . Substituting 
(C.3) into (C.2) and varying the spinor field £ we get 



5S(0 = 2 [ f>/( 

J .7=1 ^ 



M d )Bj &b d a e-eB7a b d a se _ s^^r,,. . , r + t . 



(^oid£ d ) 1/2 



where we wrote down explicitly the terms with <5£ and incorporated the terms with <5£ into the 
"c.c." (complex conjugate term). Integrating the term with daS^ by parts and taking out the 
common factor 5£ a we rewrite the above formula as 



6S(0 



/16 
(<5f)E 
,/=i 



A, 



iWjBJabd a e W]a 0M e 



+ id 



+ c.c. 



Comparing this formula with formula (7.1) we conclude that the spinor field F appearing in the 
latter is given by formula 



iW.,BJ ab d a i b 



+ id a 



wjBj ab e 

(C^«c4 d ) 1/2 



(C.l) 



where the Wj are, in turn, given by formula (C.3). 
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If we now substitute the plane wave (8.1) into formula (C.3) we get 

2£ a BJ abPaC b 



Wj = 



(C^ocdC d ) 1/2 



(C.5) 



Note that the above W.j are constant (do not depend on x), which simplifies the next step: 
substituting (8.1) into (C.4) and dividing through by the common factor e~ lpx we get 



2^> 



2WjBJ M p a ( b Wjao M C b 



(C6) 



The remarkable feature of formula (C.6) is that its RHS is constant, i.e. it does not depend on x. 
This completes the proof of Lemma 8.1. 

Let us now substitute the plane wave (8.1) directly into our Lagrangian density (C.2). Our 
Lagrangian density takes the form 



L(C; p) = J2 A -' w - 



(C.7) 



where the Wj are given by formula (C.5). The Lagrangian density (C.7) does not depend on x. 
The dynamical variable in this Lagrangian density is the constant 2-component complex spinor f , 
whereas the relativistic 4-momentum p plays the role of a parameter. Varying the spinor Q we get 



2C^0£dC d 



Wj)Wj + c.c. 



Comparing this formula with formula (8.16) we conclude that the constant spinor G appearing in 
the latter is given by formula 



G & = 2 J2 A J 



(C.8) 



It remains to observe that the right-hand sides of formulae (C.6) and (C.8) are the same. This 
completes the proof of Lemma 8.2. 



